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' In the harmonic oscillator representation, the Schrodinger equation has a form of a set of 

infinite number of algebraical equations which are labeled by the radial quantum number n. 
. It is shown that at n 3> 1 these equations are significantly simplified and become an alge- 

^ ' braical analogue of the original differential Schrodinger equation. This result is generalized 

to the case of multi-channel systems being studied in the framework of the resonating group 
, method. 

o 

^ ' §1. Introduction 

I The Schrodinger equation for a particle moving in a spherical-symmetry field 

U{r), 

{T + U{r)}^{r) = EW{r), (1) 



X 



can be transformed to the form of an infinite set of algebraical equations for the 



, Fourier coefficients of the harmonic oscillator expansion of the wave function ^'(r); 

n+l 00 

<nl\T\fil>Cni + J2<^^\Uir)\^^>Cni-ECni = 0, n, / = 0, 1, 2, ... (2) 

n=n— 1 n=0 



'Ir{r)=J2Riir)Yi^{f2), Ri{r) = J2 Cni^ni{r), (3) 



Im n=Q 



= (-')'V r(f+:!3/2) '-'^"'^^'''^' -pf--'/^'- 

where, as usually, / is the orbital momentum, n is the number of the radial quanta, 
and Lj^^^^(r^) are the Laguerre polynomials. The oscillator length ro is set to unity. 

The set of this form are obtained after the realization of the multi-channel res- 
onating group method and the hyperspherical function method in the harmonic 
oscillator basis. 

Note that even in the case of the local potential U{r) the set of its matrix 
elements 

{< nl\U{r)\fLl >,0 < n,l < 00} 



typeset using pTpJ^sty <ver.0.7> 
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is actually a discrete analogue of a non-local potential, so that the investigation of 
the infinite set (^) is complicated. In Refs.[l-3], it was shown how the set (^) can 
be completed by using the quasi-classical approximation for the Fourier coefficients, 
where one can restrict himself with a finite number of equations. However, if one 
deals with slowly decreasing potentials with power-like behavior at the infinity, the 
problem of the asymptotic Fourier coefficients remains, and for solving this problem, 
it is necessary to carry out an additional investigation. 

In this paper, the asymptotic equation in the harmonic oscillator representation 
is discussed and the asymptotic form of the potential is derived. With the use of the 
results, we can obtain a knowledge of the asymptotic behavior of a potential from its 
matrix elements. Although usually the potential itself is known, it is very important 
to apply the present method to many-body systems described with the use of a two- 
body potential. We can know the asymptotic behavior of the hyperradial potential 
even if we do not know its explicit form^. Information about the asymptotic wave 
function including the influence of long-range properties of a potential becomes very 
important in the complex scaling calculations^. 



§2. Derivation of asymptotic equations 



Let us recall how the problem of the asymptotic wave function ^{r) has been 
solved. In Eq.(|l|) with a limitation to large r, where the behavior of the potential 
U{r) becomes simpler because, in its expansion in the powers of r, the only main 
term remains, and then the precise solution of a rather simple asymptotical equation 
is found. 

In the case of the set (^), the number of radial quanta n is well analogous to the 
radius r. Therefore, we shall try to reduce those equations in (^) which correspond 
to n 1. It is known^ that at n ^ 1, the discrete variable ^J An + 21 + ?> can be 
replaced by a continuous one, r. Then, the first summation in Eq.(2), 

< nl\T\n - 1, / > C„„i,z+ < nl\T\nl > Cni+ < nl\T\n + 1,1 > Cn+i,i 

appears to be an algebraical analogue of the differential kinetic energy operator T 
multiplicated by l/\/r, i.e. 

^ \ dr dr J ^/r ' 

acting onto the wave function 



Ri = RiiVAn + 2l + 3) 



V2(4n + 2/ + 3)V4' 
Therefore, we are mainly interested in the summation 



oo 

J2 < nl\U{r)\nl > C^i. 
n=0 



(5) 
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In this summations, the terms have the following explicit form; 

/ ^nl{r)U{r)^nl{rydr / Ri{f)^nl{f)f'^df. (6) 
JO JO 

The basis of the radial functions ^nii^) satisfies the condition 

oo 

Y,'^hl{r)r<Pnl{f)r = S{r -f). (7) 

n=0 

Having changed the order of operations in (5-6) and performed summation, we have 

y2 < nl\U{r)\hl > Cfii = / <Pnl{r)U{r)Ri{f)5{r - f)rdrfdf = 
n=o -^0 

roo 

= / ^rd(r)U{r)Ri(rydr. (8) 
JO 

The basis functions 'l?ni{f) have a remarkable asymptotic behavior. That is, if 
n ^ 1, then 

where r„ = \/4n + 2/ + 3 is a turning point in the field of the harmonic oscillator 
for the particle moving with the energy 2n + / + 3/2. As a result, the last integration 
in (P) can be performed at n ^ 1, and the formula ^ is significantly simplified, 



< nl\U{r)\fil > Cni ~ U{^/An + 2l + ?,)^/2^Rl{rn) ^ [/(V4n + 2/ + 3)C7„z. 

n=0 

(10) 

Thus, in the limit of large n, the following asymptotic form of the set (^) is valid; 

n+l 

^ < nl\T\nl > Cni + UWAn + 21 + 3)C„i - ECni = 0. (11) 

h=n—l 

Bearing in mind the remark concerning the first sum in (11), we come to a conclusion 
that Eq.(ll) does not differ from the asymptotic form of the Schrodinger equation 
(1) at large r, i.e., from (11) we can turn back to (1) with a local potential U{r) but 
only if r ^ 1. 

The set ( pT] ) is a three-term recurrent relation in Cn+i^u Cn,i, and Cn-i,i- If the 
limit expression for Ri{r) at r » 1 is known, then Cn,i may be calculated due to the 
equality 

Cn,l = V2^Ri{rn). (12) 

After that, the recurrent relation (11) allows us to go backward to that values of n 
at which it remains valid. 

Having chosen this way of completing the set (^) , in the result we will have only 
small number of equations to be solved exactly. 
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§3. On the convergence to the asymptotical limit 
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Fig. 1. Matrix elements < nl\U{r)\nl > of the po- 
tential (14) with 7 = 0.25. 



Now we turn to the question; 
for which values of n the asymp- 
totic formula @) is valid? Actu- 
ally, this formula means that, at 
large n, the set of matrix elements 
< nl\U{r)\hl > is a discrete ana- 
logue of the 5-function with the nor- 
malization factor U{\/4:n + 21 + 3). 
The matrix elements themselves do 
not equal to zero only in the vicin- 
ity of the diagonal h = n, and they 
oscillate when the value of n — n is 
changed (their signs depend on the 
parity of n) , with the oscillation am- 
plitude decreasing to zero. If the 
values of U{rn) are not known a pri- 
ori, they can be found by performing 
the summation of < nl\U{r)\hl >: 

oo 

< nl\U{r)\nl > 

n=0 



= C/(V4n-h2/ + 3). (13) 

In the summation over n, the upper 
limit must be not less than 2n. Then 
it is enough to take n > 20 in order 
to make the sum converge. 

Figure 1 illustrates the be- 
havior of the matrix elements < 
nl\U{r)\nl > at different n for the 
case 

^(r) = i^^^PizZ!), (14) 
r 

while Fig. 2 (left) shows the values 
of the sums (13) for this potential. 
Convergence of the sums are ob- 
served at n = 20 and even a little 
less. 

Another example is a potential 



U{r) = — (^1 — — exp(— 7r^) 



(15) 
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which decreases as 1 /r'^ and by this reason it is characteristic of the hyperspherical 
function methodu\ The behavior of the sums (13) of the matrix elements of this 
potential is shown in Fig. 2 (right). Again, the sums are rapidly converge to the limit 
value. 

All said above remains valid for the potentials which are not singular in the 
origin. 
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Fig. 2. Asymptotics and the sums (13) of the potentials (14) (left) and (15) (right) with 7 — 0.25. 



§4. Generalization 



For a many-channel system being considered within the frame of the algebraic 
version of the resonating group method, when the reactions with light nuclei are 
investigated, the harmonic oscillator basis states \a,nl > have an additional index 
a for distinguishing between the channels. The Fourier coefficients accept the same 
index. If the basis states of different channels are coupled only by the operator of 
the potential energy, then instead of Eq.(ll) we get 

n+l p 

< «, nl\T\a, nl > C?; + J] U^p{V4n + 21 + 3)Cf; - EC^i = 0, (16) 

h=n-l 13=1 

a = 1,2, . . . ,p, 

where p is a number of channels. In other words, in this case there appear asymptotic 
matrix elements of the potential energy of the interaction between channels, and they 
contribute to the asymptotics of the Fourier coefficients C";. 

A practically inportant point (especially, in the low-energy region) is consider- 
ation of the Coulumb repulsion of the charged particles (e.g. protons) which leads 
to the polarization of clusters. Then, the matrix elements Ua/s decrease with n 
increases, 
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and it is necessary to trace over the influence of the Coulumb interaction between 
channels starting at very large n. 

The algebraic version of the hyper spherical function method employs the basis 
of the states \a,nK > where K is the hypermomentum and n includes the remain- 
ing quantum numbers. The asymptotical equations for the hyperspherical function 
methods do not differ in principle from Eqs.(lT6|). However, K appears instead of I, 
and \/4n + 21 + 3 must be replaced by y^An + 2K + 3{N - 1) where 3(iV - 1) is the 
dimension of the space in which the hyperharmonics are defined. 

§5. Conclusion 

In the harmonic oscillator representation, after the matrix elements of the Hamil- 
tonian < nl\T+U\hl > (n is the radial quantum number, I is the orbital momentum) 
are found, the Schrodinger equation has a form of a set of algebraical equations in the 
unknown Fourier coefficients Cni- In the limit of great n {n > 20), these equations 
are reduced and become a discrete analogue of the Schrodinger equation at large 
values of the radius r. Therefore, there is a one-to-one correspondence between the 
asymptotic radial function Ri{r) at large r and the asymptotic Fourier coefficients 
Cni at large n. This result can be generalized to the case of many-channel systems. 
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